The soft superconducting gap in semiconductor Majorana nanowires 
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We theoretically consider the mysterious topic of the soft gap in the tunneling conductance of the proximity 
induced superconductivity in a semiconductor-superconductor hybrid structure, where the observation of a zero 
bias conductance peak has created considerable excitement because it could be connected with the existence 
of the elusive zero-energy Majorana mode. The observed experimental superconducting tunneling gap in the 
semiconductor nanowire looks v-shaped with considerable subgap conductance even at very low temperatures 
in sharp contrast to the theoretically expected hard BCS gap with exponentially suppressed subgap conductance. 
We systematically study, by solving the appropriate BdG equations both numerically and analytically, a number 
of physical mechanisms (e.g. magnetic and non-magnetic disorder, finite temperature, dissipative Cooper pair 
breaking, interface fluctuations), which could, in principle, lead to a soft gap, finding that only the interface 
fluctuation effect is a quantitatively and qualitatively viable mechanism that is consistent with the experimental 
observations. Our work indicates that improving the quality of the superconductor-semiconductor interface, so 
that the proximity tunneling amplitude is uniform along the semiconductor nanowire, would go a long way in 
enhancing the gap in the hybrid structures being used for studying the Majorana mode. 
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Introduction. The pursuit of exotic topological phases of 
matter has become an exciting topic of research in physics |T). 
In particular, topological superconductors (SCs) supporting 
zero-energy Majorana bound states (MBS) l2T 4T"3ll have re- 
ceived increasing attention both for its intrinsic interest for 
fundamental physics as well as for its potential uses in topo- 
logical quantum computation j2l [Pfl - TTrjl . In recent proposals 
to realize topological SCs in solid state systems, an effective 
p-wave SC is induced in a semiconductor by the combined 
effects of spin-orbit coupling (SOC), Zeeman splitting of the 
energy bands and proximity induced s-wave SC (SI [TTl HHJ . 
For a semiconductor nanowire (NW), it was predicted that the 
presence of MBS at its ends could be experimentally detected 
in the differential tunneling conductance G ( V) = dl/dV at the 
interface with a normal contact via the emergence of a zero- 
bias peak (ZBP) of height 2e 2 /h (at zero temperature) |[T9l - 
[23). Interestingly, recent experiments implementing the above 
NW proposal using InSb or InAs NWs in contact with a SC 
(mostly Nb) have shown a suggestive ZBP in the G (V) spec- 
tra that is consistent with the theoretical predictions of the 
MBS scenario [24-26]. However, despite the enormous ex- 
citement generated by these results there are several discrep- 
ancies between theory and experiment that have attracted re- 
cent attention, e.g. the suppressed value of the ZBP compared 
with its predicted quantized Majorana value of 2e 2 /h 11271 . 
the non-closing of the superconducting gap at the topological 
phase transition in the finite magnetic field 11281 . and the short 
wire length leading to possible Majorana overlap between the 
ends ||29ll . In addition, it has been pointed out that other mech- 
anisms, notably disorder Il30ti33l and Kondo effect [34], could 
possibly lead to a ZBP in the experiments. 

Curiously, an ubiquitous (and rather striking) feature of the 
experiment involving proximity-induced superconductivity in 
the semiconductor itself has remained ignored in the literature 
despite a great deal of activity in the subject: the supercon- 



ducting gap in the NW, as reflected in the measured G(V), 
is extremely "soft" in both the high-field topological phase 
(where the ZBP exists) and in the zero-field or the low-field 
trivial phase where there is no ZBP. In fact, the soft gap fea- 
ture, which is clearly a property of the semiconductor-SC hy- 
brid structure quite independent of the MBS physics, is ex- 
tremely prominent in the data (even in the absence of the ap- 
plied Zeeman field) with the subgap conductance being typ- 
ically only a factor of 2-3 lower than the above-gap conduc- 
tance, implying the existence of rather large amount of subgap 
states whose origin remains unclear. Amazingly, this soft gap 
issue has not been discussed at all in the theoretical literature, 
in spite of the soft gap being present in all of the experimental 
data from the groups so far appearing in the subject, and all 
the theoretical calculations invariably manifest the usual BCS 
hard gap (with the Majorana ZBP superposed on this hard gap 
above the magnetic -field-induced topological quantum phase 
transition). We believe that without a thorough understanding 
of this ubiquitous soft gap feature of the experiment the whole 
subject remains incomplete and even highly suspect. 

In this Letter we provide the first theoretical attempt to ex- 
plain the observed soft gap and perform a systematic theoreti- 
cal investigation of different possible mechanisms giving rise 
to subgap tunneling conductance using a theoretical model 
similar to the experimental configurations in Refs. l24] - l2"6l In 
particular, we focus on the effects due to: (a) non-magnetic 
and (b) magnetic disorder in the NW; (c) temperature; (d) dis- 
sipative quasiparticle broadening arising due to various pair- 
breaking mechanisms such as poisoning, coupling to other de- 
grees of freedom (e.g. phonons) or due to electron-electron 
interactions; and (e) inhomogeneities at the SC-NW interface 
due to imperfections (e.g. roughness and barrier fluctuations) 
that may arise during device fabrication. Since the soft gap 
occurs universally in the experiment at all parameter values, 
we consider only the non-topological zero-magnetic field situ- 
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ation in this paper because this is where the gap should be the 
largest and the hardest. Our results for the topological phase 
with finite magnetic field (not shown in this paper) are very 
similar qualitatively. Our numerical results point to a stronger 
effect of inhomogeneities at the semiconductor NW-SC inter- 
face (due to the exponential dependence on the barrier param- 
eters), as compared to disorder, temperature or pair-breaking 
broadening effects, in producing the soft gap as described in 
the rest of this paper. We substantiate our numerical calcu- 
lations with an analytical approach using Abrikosov-Gor'kov 
theory [35 1 within a minimal model for a semiconductor NW 
under the effect of a fluctuating s-wave pairing potential in- 
duced at the interface. 

Theoretical model. We consider the Hamiltonian describ- 
ing a narrow one-dimensional semiconductor NW of length 
L x along the x-axis subject to SOC, Zeeman field along its 
axis, and proximity-induced s-wave pairing, 

H K = J <fx|cj"(x) - ^ - /i(x) + ia R o- y d x - B z cr x 

-b(x)-a\ c s ,(x) + A(x)(cl(x)c|(x) + H.c.)l. (1) 

- ss' I 

Here, c\(x) creates an electron with spin projection s =T, 1 at 
position x in the NW, m* is the effective mass of the electrons 
in the semiconductor, a« is the Rashba velocity related to the 
SOC energy E so = m* e al/2, B z is the Zeeman energy due to 
an external magnetic field along, <x = (tr x , tr y , cr z ) is the vec- 
tor of Pauli matrices, and A(x) is the (fluctuating) proximity- 
induced s-wave SC pairing. The sums over spin indices are 
implied. Static non-magnetic disorder in the semiconduc- 
tor NW is included through a fluctuating chemical potential 
//(x) = fj.Q+5/u(x) around the average value /jq. Static magnetic 
disorder may be present in the sample due to contamination 
with magnetic atoms or due to the presence of regions in the 
NW acting as quantum dots with an odd number of electrons. 
Assuming that the experimental temperature T is higher than 
the Kondo temperature T K associated with the magnetic impu- 
rities, we can neglect the quantum dynamics of the impurity- 
spins and model its effect as a randomly oriented inhomoge- 
neous magnetic field b(x) |36|. Finally, we also introduce spa- 
tial (and associated potential) fluctuations in the barrier sep- 
arating the NW and the SC. These fluctuations give rise to 
spatial fluctuations in the induced s-wave SC pairing in the 
NW due to proximity effect. The induced pairing is very well 
approximated by the relation A(x) = y(x)Asc/ (y(x) + Asc), 
where Asc is the uniform parent pairing in the bulk SC, 
and y(x) = pQt]_(x) is the transparency of the NW-SC in- 
terface 1371 . Here, po is the local density of states of elec- 
trons in the NW at the Fermi energy in the normal phase, and 
fj_ (x) is the tunneling matrix element connecting the SC and 
NW. The effect of the interface inhomogeneities is modeled 
through t x (x) = t\e~ KSd ^ x \ where 6d(x) denotes the fluctua- 
tion in the width of the NW-SC barrier and k is a phenomeno- 
logical constant with units of inverse length that parametrizes 
the energy barrier of the NW-SC interface. Such a functional 



form is expected due to fluctuations in the overlap of evanes- 
cent wavefunctions. For y(x) <K Asc, the induced SC pair- 
ing is controlled by the interface transparency A(x) w y(x) 
113711 , and thus can be modeled as A(x) = Aoe~ 2s ^ x \ where 
the dimensionless parameter 6B (x) = Kdd (x) characterizes the 
roughness of the interface, and Ao is the induced SC pairing 
in the absence of the interface inhomogeneity. Note that our 
model for interface fluctuations is generic and only incorpo- 
rates the invariable presence of potential fluctuations at the 
interface separating the SC metal and the NW. A good inter- 
face by definition will have very little fluctuations whereas a 
poor interface will have large fluctuations. 

We approximate ([TJ by a tight-binding Hamiltonian defined 
on an A^-site lattice, 

an i 
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where a = au/2a, and a is the lattice constant. The estimated 
parameters for the InSb NW are L x = 2pm, m* = 0.015m e , 
£.50 = m* e a 2 R l2 = 50/ieV, Ao = 250^eV, and temperature 
T = 70mK a 0.024A E3). For our tight-binding simula- 
tions we use a one-band model with = 500, t - 616fieV, 
and a = 0.07f. The different disorder mechanisms are taken 
into account by introducing Gaussian-distributed random vari- 
ables d>yU,, bj - (b x ,b y n bfj, and 5/3, with zero mean and vari- 
ance given by {SmSfijJ = W 2 S U , (tfb*) = Wl8ij8 pq , and 

(dB/dB/j = Wpdij, respectively. To model the interface inho- 
mogeneity, we coarse-grain the interface in patches of length 
5a and assume that 5B, is uniform within each patch, but varies 
randomly from patch to patch with a standard deviation of Wp. 
Note that assuming a Gaussian distribution in SBi results in a 
different probability distribution function for A ( , 
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from which the mean (A,) = dA, P(A,)A ; and variance 
W 2 h = §™ dAj P (A,) (A,- - (A,)) 2 are calculated. For example 
a typical variation of 10% in the width of the barrier, i.e. Wp = 
0.1, causes a variation W&/ (A,-) ~ 0.2 in the s-wave pairing. 

The relevant experimental quantity is the tunneling differ- 
ential conductance G (V) at an end of the NW, which is related 
to the local density of states 11241 1251. We calculate G(V) us- 
ing the tunneling formalism by coupling the NW to a contact 
lead 12T1 I22 38 1 . The Hamiltonian of the combined system is 
H - H K + Hi + H t , where Hi = 2*j £kd\ s dk s is the Hamil- 
tonian describing the lead and H t = ti Y*ks c hs Cls + ^- c - ^ s ^ e 
tunneling Hamiltonian coupling site / — 1 of the NW to the 
lead via a tunneling matrix element t L . The tunneling conduc- 
tance at the site i = 1 is then 

dajp'l (cj) f (co - eV) , (4) 

OO 
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FIG. 1. (color online) Differential conductance for electron tunneling 
into an end of the semiconductor nanowire for zero external magnetic 
field. Various possible scenarios for subgap states are considered: 
(a) static disorder, (b) magnetic disorder, (c) thermal pair-breaking 
scenario and (d) quasiparticle broadening scenario. See text for more 
details. 

where f(x) = (e x ^ T + is the Fermi distribution function, 
pi is the density of states at the Fermi energy in the lead, and 
V is the voltage at which the lead is biased with respect to /jq. 
Here, p™ (cS) is the local density of states in the NW (including 
both spin projections) at site i = 1 in the presence of the lead, 
which we calculate as p!" (a>) = — ^Im (w). Here g™. (oS) is 
the retarded Green's function of the NW in real-space repre- 
sentation, which in the limit 1 1 — » becomes 

(0)* (0) (0)* (0) 

Zu\ u . v. v ■ 
is,n js,n is,n js,n 
— + T F( o) T ■ — ' (5) 

with and [uf^ , vf? } being, respectively, the eigenvalues 
and eigenvectors resulting from the diagonalization of the 
BdG Hamiltonian corresponding to Eq. (|2]). To include the 
presence of the lead, we solved the equation of motion for 
g™j (z) in the presence of the tunneling term H t ||39l . The term 
jm is the self-energy, which in the limit ti — > becomes 
y L ,n = -in PL tl^ufJ + \vfj). 

Results. We now present the numerical results for G ( V). A 
mean chemical potential of fio = -338//eV is used, and the 
temperature is set to T — 70mK ll24ll unless otherwise stated. 
In Fig.[TJa) we consider the effect of static disorder on G(V). 
We take W^ — (blue curve) to = 0.8Ao (red curve) in 
equal steps of 0.2Ao. The plots are offset in steps of 0.1 for 
clarity. As expected from Anderson's theorem 061 l40l . our 
results show that the subgap density of states is not affected 
by the presence of static non-magnetic disorder, thus render- 
ing this an unlikely mechanism for the observed subgap con- 
ductance. We note as an aside that in our numerical results for 
the topological phase, which are not shown here, the effect of 
non-magnetic disorder is stronger than in the zero magnetic 



field non-topological phase since Anderson's theorem does 
not apply in the topological phase. In fact, the non-magnetic 
disorder in the topological phase behaves very similar to the 
magnetic disorder in the non-topological phase discussed be- 
low. 

The effect of magnetic disorder is shown in Fig. |TJb). We 
have taken W b = 0, 0.27 Ao, 0.54A , 0.68A and 0.8 1A (blue 
to red curves). In this case, we find a substantial modification 
in the subgap conductance. In particular, a soft superconduct- 
ing gap, similar to the one observed in Ref. 24, is obtained 
for Wb = 0.8 lAo (red curve). Although we cannot exclude 
the presence of magnetic disorder, a large amount is needed 
to explain the data. We obtain an estimate A0T1, ~ 1 where 
Tj, = 2f 2 (l - (///2f) 2 )/3vpW 2 is estimated from our tight bind- 
ing parameters. Such a large amount of magnetic disorder is 
unlikely to be present in the NW used in the experiments. 

The thermal pair-breaking effect is considered in Fig.[TJc). 
Note that in our non-interacting theoretical description, ther- 
mal effects can only arise from broadening of the Fermi dis- 
tribution function in Eq. Q. We vary the temperature from 
T - 0.027 A () (blue curve) to T — 0.35Ao (red curve) in equal 
steps of 0.054A (the conversion is 0.027 A ~ 78mK). Al- 
though a considerable amount of thermally-induced subgap 
conductance is obtained for T = 0.35Ao (red curve), this 
value is much larger than the reported experimental temper- 
ature r exp = 70 mK, and cannot by itself explain the exper- 
imental features. We note that the blue curve corresponds to 
T = 78mK > r exp , for which there is no appreciable subgap 
conductance. 

In Fig.[TJd), we consider the effect of a finite quasiparticle 
broadening due to dissipative processes. While it is exper- 
imentally unclear what physical mechanism would give rise 
to this broadening, it can in principle arise due to coupling of 
electrons in the NW to a source of dissipation, e.g. presence of 
(unconsidered) normal contacts, quasiparticle poisoning due 
to tunneling of normal electrons into the NW, scattering with 
phonons and/or other electrons. These pair-breaking mecha- 
nisms can be phenomenologically modeled by a shift in the 
frequency a> — > to + iy^ in Eq. ([5]), where is the quasi- 
particle broadening. Quasiparticle lifetime effects were con- 
sidered in a similar way in the context of BCS superconduc- 
tors by introducing a phenomenologically broadened density 
of states p(to,y N ) = Re[(\to\ + iy N )/[(\oj\ + iy N ) 2 + A 2 ] l/2 ] ED- 
In Fig. |TJd) we vary y^ from y^ = 0.027Aq (blue curve) to 
yN - 0.35Ao (red curve) in equal steps of 0.054Ao. We see 
that even for the largest values of y^ (i.e. y^ ~ 0.35Ao corre- 
sponding to the red curve), a remnant of the hard SC gap is still 
present. Therefore, this effect alone is incapable of explaining 
the substantial gap softening observed in the experiments. 

While all of the above-mentioned mechanisms are likely to 
be present to some extent in a realistic setup, our results in- 
dicate that it is unlikely that they can individually explain the 
experimentally observed soft gap. Moreover, even after com- 
bining all the effects of non-magnetic and magnetic disorder, 
quasiparticle decay rate of order 0.1 Ao, and temperature of 
70mK, we found that obtaining a soft gap that qualitatively 



FIG. 2. (color online) Differential tunneling conductance for zero 
external magnetic field in the presence of SC-NW interface disorder 
and quasiparticle broadening. Two different interface inhomogeneity 
strengths are taken: (a) Wa = 0.4 and (b) Wp = 0.8. In each plot, we 
consider quasiparticle decay rates from y N = 0. 1 1A (blue curve) to 
jn = 0.32Ao (red curve) in equal steps of 0.054Ao. See text for more 
details. 

agrees with experiments requires magnetic disorder strength 
of AoTfo ~ 0(1), which seems to be unrealistic. This leads 
us finally to the effect of inhomogeneities at the NW-SC in- 
terface (see Fig. pi. More generically, pair-breaking mecha- 
nisms at the interface have also been put forward to explain 
the transport properties in other semiconductor-SC systems 
that show similar subgap conductance |42|. We now argue 
that a reasonable amount of interface inhomogeneity together 
with quasiparticle broadening gives a soft gap that is in good 
qualitative and semi-quantitative agreement with the experi- 
mental findings, thus rendering the combination of these two 
effects as the most likely candidate for the soft gap. In Fig. [2] 
we take (a) Wp = 0.4 and (b) Wp = 0.8. In each plot, the 
five curves correspond to y^ varying from y^ = 0.1 lAn (blue 
curve) to y^ = 0.32Ao (red curve), in equal steps of 0.054Aq. 
We observe a large amount of subgap contributions, with a 
noticeable "v-shaped" tunneling conductance around V = 0, 
which is consistent with experimental observations l24l l25l . 
In Fig. |2|b), we see that y^ ~ 0.1 Ao is sufficient to obtain 
a soft gap reminiscent of the experimental findings. We add 
that the v-shaped soft gap is obtained only in the presence 
of both the interface fluctuations and quasiparticle broaden- 
ing, and an unrealistic magnitude for either of these depairing 
mechanisms is needed to reproduce the soft gap in the absence 
of the other. 

An order-of-magnitude estimate for the dimensionless pa- 
rameter 6/3, can be obtained based on known experimental pa- 
rameters. The width of the NWs used in Ref. 24 was quoted 
as lOOnm ± lOnm. Assuming that the fluctuations in the SC- 
NW barrier width is of order the wire width fluctuations, we 
take Sdj ~ 5nm. The phenomenological barrier parameter 
k can be estimated using the interface energy barrier Uq via 
k * ^2m* e U()/h. Using an estimate for Uq based on a Nb- 
InGaAs junction [43], we take Uq « 0.2 eV. With an effective 
mass for the InSb wire, m* « 0.015ra c , we obtain 5/3j ~ 1. 
This order of magnitude estimate is consistent with the stan- 
dard deviations Wp = 0.4 and 0.8 used in this work. 

A minimal analytical model that captures the relevant ef- 
fects can be obtained from Eq. ([T| in the absence SOC, Zee- 



FIG. 3. (color online) Analytical results for p'*{a))/po, the averaged 
local density of states obtained from the self-consistent Abrikosov- 
Gor'kov for various values of A t a , and (a) y N = and (b) y N = 
0.1 lAo. See text for more details. 

man effect and other types of disorder, and assuming the SC 
pairing itself to be a Gaussian variable A(x) = Aq + SA(x) with 
variance (SA(x)dA(x')) = W^6 (x - x'). This model allows 
one to obtain a useful insight into the effects of a fluctuating 
SC pairing on G(V). We use the theoretical framework of the 
self-consistent Abrikosov-Gor'kov theory f35l [36l to obtain 
the averaged electron Green's function gfj(oj). The calcula- 
tions are shown in detail in the supplementary material Il44l . 
Despite the mathematical similarity of the formalism to the 
(more usual) case of scattering induced by magnetic impuri- 
ties in s-wave SCs, here we are only considering SC pairing 
fluctuations as the pair-breaking mechanism. In Fig.[3|a), we 
show the results for pj V (w)/po, the averaged local density of 
states at the end of the NW, which is the main quantity deter- 
mining G(V) at T — [cf. Eq. Q]. In each plot, the black to 
purple curves correspond to (AoTa) -1 = to 1.5 in equal steps 
of 0.25. Here, r A ' = nW^po is the scattering rate induced by 
SC pairing fluctuations. Interestingly, the theory allows us to 
obtain an analytical expression for the quasiparticle gap spec- 
trum in p™(u))/p : £ gap = Ao[l - (A t a )- 2 / 3 ] 3 / 2 El. For 
Aota < 1, the quasiparticle gap vanishes (brown curve). To 
make contact with our numerical results in Fig. [2] in Fig.|3jb) 
we consider a finite y N = 0.1 IAq, for the same values of Aqt a 
as in Fig.[3ja). The quasiparticle decay rate y^ has the effect 
of broadening the sharp edge features present in the tunnel- 
ing spectrum when y^ - 0. Again, we see that fluctuations 
in the induced SC pairing together with quasiparticle broad- 
ening gives the characteristic v-shaped tunneling spectrum in 
the subgap regime (e.g. cyan and green curves). Despite the 
minimal nature of our analytical model, the theory qualita- 
tively explains the results of our numerical simulations. 

We note that pairing fluctuations in the parent SC may also 
play a role here since they will also induce pairing fluctuations 
inside the nanowire l45l . However, given the universality of 
the soft gap behavior in semiconductor-SC hybrid structures, 
which appears independently of the material being used for 
the parent SC, and that we are in the weak-coupling regime 
(i.e. y(x) <K A$c), we believe that the soft gap behavior is 
mainly caused by the interface fluctuations. 

To summarize, we have studied the effect of different pair- 
breaking mechanisms likely arising in semiconductor-SC Ma- 



5 



jorana NWs, and systematically analyzed their influence on 
the subgap tunneling conductance in order to explain the ex- 
perimentally observed soft gap behavior. While we cannot 
completely rule out some of these mechanisms (i.e. mag- 
netic scattering, thermal and dissipative broadening), quan- 
titative considerations point to the interface fluctuations at 
the semiconductor-SC contact leading to inhomogeneous pair- 
ing amplitude along the wire as the primary physical mech- 
anism causing the ubiquitous soft gap behavior. Our work 
indicates that materials improvement leading to optimized 
semiconductor-SC interfaces should considerably ameliorate 
the proximity gap in the hybrid structures. 

The authors are grateful to Meng Cheng and Claudia Ojeda- 
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from DARPA QuEST, JQI-NSF-PFC and Microsoft Q. 



[I] A. Kitaev, Usp. Fiz. Nauk 171, 131 (2001), supplement; Pro- 
ceedinds of Chernoglovka 2000, cond-mat/001440. 

[2] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and 
S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

[3] S. Das Sarma, C. Nayak, and S. Tewari, Phys. Rev. B 73, 220502 
(2006). 

[4] S. Tewari, S. D. Sarma, C. Nayak, C. Zhang, and P. Zoller, Phys. 

Rev. Lett. 98, 010506 (2007). 
[5] L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008). 
[6] C. Zhang, S. Tewari, R. M. Lutchyn, and S. Das Sarma, Phys. 

Rev. Lett. 101, 160401 (2008). 
[7] M. Sato and S. Fujimoto, Phys. Rev. B 79, 094504 (2009). 
[8] J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, Phys. 

Rev. Lett. 104, 040502 (2010). 
[9] R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. Lett. 

105, 077001 (2010). 
[10] Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 105, 

177002 (2010). 

[II] J. Alicea, Phys. Rev. B 81, 125318 (2010). 

[12] L. Jiang, T. Kitagawa, J. Alicea, A. R. Akhmerov, D. Pekker, 

G. Refael, J. I. Cirac, E. Dernier, M. D. Lukin, and P. Zoller, 

Phys. Rev. Lett. 106, 220402 (2011). 
[13] L. P. Rokhinson, X. Liu, and J. K. Furdyna, 

10.1038/nphys2429. 
[14] J. D. Sau, S. Tewari, and S. Das Sarma, Phys. Rev. A 82, 052322 

(2010). 

[15] J. Alicea, Y. Oreg, G. Refael, F von Oppen and M. P. A. Fisher, 

Nature Physics 7, 412 (201 1). 
[16] B. van Heck, A. R. Akhmerov, F. Hassler, M. Burrello, and C. 

W. J. Beenakker, New J. Phys. 14 035019 (2012). 
[17] R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. Lett. 

105, 077001 (2010). 
[18] Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 105, 

177002 (2010). 

[19] K. Sengupta, I. Zutic, H.-J. Kwon, V. M. Yakovenko, and 

S. Das Sarma, Phys. Rev. B 63, 144531 (2001). 
[20] K. T. Law, P. A. Lee, and T. K. Ng, Phys. Rev. Lett. 103, 237001 

(2009). 

[21] J. D. Sau, S. Tewari, R. M. Lutchyn, T. D. Stanescu, and 

S. Das Sarma, Phys. Rev. B 82, 214509 (2010). 
[22] K. Flensberg, Phys. Rev. B 82, 180516 (2010). 
[23] M. Wimmer, A. R. Akhmerov, J. P. Dahlhaus, and C. W. J. 



Beenakker, New J. Phys. 13, 053016 (2011). 
[24] V. Mourik, K. Zuo, S. M. Frolov, S. Plissard, E. A. Bakkers, 

and L. Kouwenhoven, Science 336, 1003 (2012). 
[25] A. Das, Y. Ronen, Y. Most, Y. Oreg, M. Heiblum, and H. Shtrik- 

man, ArXiv e-prints (2012), 1205.7073. 
[26] C. M. Marcus, private communication. 

[27] C.-H. Lin, J. D. Sau, and S. Das Sarma, ArXiv e-prints (2012), 
1204.3085. 

[28] T. D. Stanescu, R. M. Lutchyn, and S. Das Sarma, ArXiv e- 

prints (2012), 1208.4136. 
[29] T. D. Stanescu, S. Tewari, J. D. Sau, and S. Das Sarma, ArXiv 

e-prints (2012), 1206.0013. 
[30] D. Bagrets and A. Altland (2012), arXiv: 1206.0434. 
[31] J. Liu, A. C. Potter, K. Law, and P. A. Lee (2012), 

arXiv: 1206. 1276. 
[32] D. I. Pikulin, J. P. Dahlhaus, M. Wimmer, H. Schomerus, and 

C. W. J. Beenakker (2012), arXiv:1206.6687. 
[33] G. Kells, D. Meidan, and P. W. Brouwer (2012), 

arXiv: 1207.3067. 
[34] E. J. H. Lee, X. Jiang, R. Aguado, G. Katsaros, C. M. Lieber, 

and S. De Franceschi, Phys. Rev. Lett. 109, 186802 (2012). 
[35] A. A. Abrikosov and L. P. Gor'kov, Zh. Eksp. Teo. Fiz. 39, 178 1 

(1960). 

[36] A. V. Balatsky, I. Vekhter, and J.-X. Zhu, Rev. Mod. Phys. 78, 
373 (2006). 

[37] T. D. Stanescu, R. M. Lutchyn, and S. Das Sarma, Phys. Rev. B 

84,144522(2011). 
[38] Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 (1992). 
[39] G. D. Mahan, Many-Particle Physics, Physics of Solids 

and Liquids (Kluwer Academic/Plenum Publishers, New York, 

2000), 3rd ed. 

[40] P. W. Anderson, J. Phys. Chem. Solids 11, 26 (1959). 

[41] R. C. Dynes, V. Narayanamurti, and J. P. Garno, Phys. Rev. 

Lett. 41, 1509 (1978). 
[42] K. Neurohr, A. A. Golubov, T. Klocke, J. Kaufmann, 

T. Schapers, J. Appenzeller, D. Uhlisch, A. V. Ustinov, 

M. Hollfelder, H. Liith, et al„ Phys. Rev. B 54, 17018 (1996). 
[43] A. Kastalsky, A. W. Kleinsasser, L. H. Greene, R. Bhat, F. P. 

Milliken, and J. P. Harbison, Phys. Rev. Lett. 67, 3026 (1991). 
[44] See Supplemental Material for detailed calculations. 
[45] E. F. C. Driessen, P. C. J. J. Coumou, R. R. Tromp, P. J. 

de Visser, and T. M. Klapwijk, Phys. Rev. Lett. 109, 107003 

(2012). 



6 



Appendix: Supplemental Materials for "The soft 
superconducting gap in semiconductor Majorana nanowires'' 



For notational convenience we rewrite the fluctuation- 
averaged Green function as 



In this supplementary document we present the details of 
the Abrikosov-Gor'kov theory [ 1 1 adapted to handle pairing 
potential fluctuations. 

Let us consider the Dyson's equation for the Green function 
in the presence of fluctuating pairing potential: 

Gkk' = Gf Skkf + / GfWkk'G^ 



P 



'k! 
MO) 



+ ]>]Gf%,GfV pk ,G^ (A.1) 
p 

+ Y,GTv kp GfV pp ,G^V p , k ,Gf + ... (A.2) 



/ , 'k 

PP' 



where G^ 0) = (zto - &t 3 - A0T1) _1 is the Green function 
for the superconducting wire with uniform Ao, and V kp - 
(6A) k+p Ti is the fluctuation in pairing potential, r, are the 
Pauli matrices acting on the particle-hole space. If we con- 
sider a Gaussian white noise for A : 



(6A k ) = 
(6A k 6A t ) = W 2 A 6 k ,- t 



(A3) 
(A.4) 



then the fluctuation-averaged Green function can be evaluated 
with cumulant expansion: 

G k 5 kt = Gf W + Yj Gf TiG^nG® (6A kp 6A pt ) 



(0) 



/ 1 k 

PP' 



Gk-G k + W A G k T] 



(5A kp 6A pp ,5A p , p „6A p „ k ,) + . 



(A.5) 



t r (0) 
T iG k 



+ WlGf Tl 



\ P ' 

v p 



2>r 



V P 

Y.Gfr^Gf 

Vpp' 



Zwo), r to) r jm 
G p T \G p , TiG k+p ,_ p 



TlG 



(0) 



'k 



TlG 



CO) 



(A.6) 



Note that this expression is formally exact, provided we in- 
clude all higher order processes. We now consider the follow- 
ing approximation for the self-energy: I. k - W A J ^T\G P T\. 
This corresponds to the "non-crossing approximation" where 

terms like Y,„ n , G\?T\G^}r\GV > , are ignored. This is a 
t-tpp p 1 p 1 1 k + P f - P & 

good approximation if k F v F r A » 1, where r A ' = nW^po is 

the scattering rate and po is the density of states at the Fermi 

point for the wire with uniform A. 



G k = (ZT - &T 3 - A T! - E) 1 

= (zT -|jfcT 3 - An) -1 



(A.7) 



where z = Z — Eo > si = fit + ^3 an d A = Ao + Si . Here the self- 
energy is split into its components: E = E0T0 + E3T3 + E1T1 

The self-energy can be readily evaluated as 
' dp z + #*t 3 + Ati 



In z 
-1 z + Ati 



2 -e P 



A 2 



2t a ^2 



(A.8) 



Substituting this expression into Eq. (A.7i yields the renor- 
malization of the energy and pairing potential: 



1 



z + 



A = An 



2t a Jfc 
1 



2r A ^2 _ ~ z 2 



(A.9) 
(A. 10) 



The density of states in the presence of pairing potential 
fluctuation is obtained by 



P(e) = -^Im j 

1, 



2n 



TvG k (z -> e + iO) 



= po (A t a ) Imw (A. 11) 

where we have defined u = z/A and the relation 



An 



is used in the last step. The gap in the quasi- 



particle spectrum is the smallest a> at which u acquires an 
imaginary part. For the function z(u) for real u, the maximum 
value of z, beyond which u becomes complex, is Q: 

£ gap = 0(A q t a - 1) A (l - (A r A r 2/3 ) 3/2 (A.12) 

To compute the local density of states near the end of a 
semi-infinite wire, we add an impurity at r — of an infinite 
wire: Vi mp = U6{r)T3. Taking the limit U — > 00 will cut the 
wire into two halves. Consider now the Dyson's equation for 
the fluctuation-averaged Green function in the presence of the 
impurity: 

G irap (r, r) = G(r, r) + G (0, r) C/r 3 G imp (0, r) 
G imp (0, r) = G(0, r) + G (0, 0) t/r 3 G lmp (0, r) 

= [l -G(0,0) f/T 3 ] _1 G(0,r) 
.-. G imp (r, r) = G (0) + G (r) r 3 [ir 1 - G (0) t,]^ G (-r) 
where G(r) = J ^e' kr G k . Hence the Green function in the 



presence of a boundary is Gb (r, r) = Gj mp (r, r) 



G(0)- 
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G (r) G (0) 1 G (-r). The local density of states at r is then 

^ImTrGb (r, r). 
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